Strange quark asymmetry in the nucleon and the NuTeV anomaly 
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The NuTeV anomaly of a non-universal value of the fundamental parameter sin 2 9w in the elec- 
troweak theory has been interpreted as an indication for new physics beyond the Standard Model. 
However, the observed quantity depends on a possible asymmetry in the momentum distributions 
of strange quarks and antiquarks in the nucleon. This asymmetry occurs naturally in a phenomeno- 
logically successful physical model for such parton distributions, which reduces the NuTeV result to 
only about two standard deviations from the Standard Model. 
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The Weinberg angle 9w is a fundamental parameter 
in the electroweak theory pj. It quantifies the 'mixing' 
of the electromagnetic and weak forces resulting in the 
photon and Z° as physical quanta mediating these in- 
teractions. Within the theory, the value of 9w should 
be universal for all possible observables. It is, therefore, 
of great interest that the value of sin 2 9\y recently ex- 
tracted from measurements of neutral and charged 
current neutrino and anti-neutrino cross-sections by the 
NuTeV collaboration is about three standard deviations 
above the value obtained from previous measurements, 
e.g. from e + e~ annihilation at LEP. A number of possible 
explanations Q for this discrepancy have been suggested 
in terms of extensions to the Standard Model, but also 
in terms of inadequate descriptions of effects within the 
Standard Model. 

An explanation of the latter kind would be a difference 
in the momentum distributions of strange (s) and anti- 
strange (s) quarks in the nucleon, since neutrinos and 
anti-neutrinos interact differently with s and 5. Although 
such ss pairs occurring as quantum fluctuations must 
have equal numbers of s and s to conserve the strangeness 
quantum number, their momentum distributions may be 
different. The experimental evidence for such an asym- 
metry is as yet inconclusive l| El 0> Qj but there are 
theoretical motivations for it In this paper we show 
that an asymmetric strange sea is a natural consequence 
of a physical model for parton momenta that reproduces 
experimental data on the proton structure, such as the 
structure function and the u — d asymmetry. In our 
model, the significance of the NuTeV anomaly is reduced 
to about two standard deviations, leaving no strong hint 
for new physics beyond the Standard Model. 

In the NuTeV experiment, the value of sin 2 9w is ex- 
tracted from neutral and charged current cross-sections 
of neutrinos and anti-neutrinos, using a procedure based 
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on the Paschos-Wolfenstein relation 0, ITo| 
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where the neutral current quark couplings g are given by 
2 V- Here, r = a(PN -> l + X)/a{vN -> i~X) ~ 1/2 
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are the primary observables of NuTeV. The small devi- 
ation from an isoscalar target by an excess of neutrons 
over protons is accounted for in the NuTeV analysis, but 
exact isospin symmetry in neutron and proton quark dis- 
tributions is assumed_, i.e. u p (x) = d n (x), d p (x) = u n (x) 
and similarly for u, d. As usual, x represents the frac- 
tional momentum of the quark in the nucleon. The re- 
sult for R~ in eq. is, furthermore, based on symme- 
try in the strange and charm quark sea distributions, i.e. 
s(x) = s(x) and c{x) — c(x). If these symmetries are 
not present in the nucleon, the NuTeV result includes an 
erroneous shift in the value of sin 2 9w Q- • 

For the part of the nucleon sea arising from gluon 
splittings g — > qq in perturbative QCD, symmetry is ex- 
pected in the distributions of quarks and antiquarks, i.e. 
q(x) — q(x). Conventional parameterisations of quark 
momentum distributions assume this symmetry also for 
the cc-distribution at the start of the perturbative QCD 
evolution. However, for these sea distributions arising 
from the non-perturbative dynamics of the bound state 
nucleon there may well be such asymmetries. We show 
that nucleon fluctuations into \AK), where the s quark 
is in the heavier A baryon and the s is in the lighter K 
meson, give a harder momentum distribution for the s 
than the s which thus affects the NuTeV analysis. 

We have previously [TllllJI presented a physical model 
giving the momentum distributions of partons in the nu- 
cleon, as illustrated in Fig. ^ More precisely, the model 
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Figure 1: Probing a valence parton in the proton and a sea 
parton in a hadronic fluctuation (letters are four- momenta). 

gives the x-shape of the parton distributions at a mo- 
mentum transfer scale Qq w 1 GeV 2 , i.e. xq(x,Ql) and 
xg(x, Qq), which provide an effective description of the 
non-perturbative dynamics of the bound state nucleon. 
Our approach gives the four-momentum k of a single 
probed parton, whereas all other partons are treated col- 
lectively as a single remnant with four-momentum r cor- 
responding to integrating out all other information in the 
nucleon wave function. It is convenient to consider the 
nucleon rest frame where there is no preferred direction 
and hence the parton momentum distribution is spheri- 
cally symmetric. The shape of the momentum distribu- 
tion for a parton of type i and mass rrii is then taken as 
Gaussian 

f i (k)=N{a i ,m i )exp{ — — J ( 3 ) 

which may be motivated as a result of the many inter- 
actions binding the parton in the nucleon. The width of 
the distribution should be of order hundred MeV from 
the Heisenberg uncertainty relation applied to the nu- 
cleon size, i.e. = 1/dpj. This Fermi motion inside the 
nucleon provides the 'primordial transverse momentum', 
which has been extracted from deep inelastic scattering 
data and found to be well described by a Gaussian dis- 
tribution of a few hundred MeV width giving phe- 
nomcnological support for this description. 

The energy component does not have the same simple 
connection to the Heisenberg uncertainty relation. It is 
assumed to have the same Gaussian fluctuation around 
the parton mass, such that partons can be off-shell at 
the soft scale of the binding interactions. This means a 
parton fluctuation life-time corresponding to the nucleon 
radius. 

The momentum fraction x of the parton is then de- 
fined as the light-cone fraction x = k + /p + . Here, four- 
momenta are expressed as p = (p + ,p-,p±) where the 
'plus' and 'minus' components are p± = E ±p z and the 
z-axis defined by the probe. The fraction x is then invari- 
ant under boosts along the z-axis and equivalent to the 
conventional momentum fraction x = k z /p z in a frame 
where p z is large ('infinite momentum' frame). 

In order to obtain a kinematically allowed final state, 
one must impose the following constraints. The scat- 
tered parton must be on-shell or have a time-like vir- 
tuality (causing final state QCD radiation), i.e. have a 



mass-squared in the range m| < j 2 < W 2 (W is the in- 
variant mass of the hadronic system). Furthermore, the 
hadron remnant r is obtained from energy-momentum 
conservation and must have a sufficient invariant mass 
to contain the remaining partons. These constraints also 
ensure that < x < 1. 

The parton distributions are obtained by integrating 
eq. © with these conditions. Using a Monte Carlo 
method this can be achieved numerically without approx- 
imations. (With further simplifying approximations it is 
possible to obtain an analytical expression for the parton 
densities H3-) The normalisation of the valence distri- 
butions is provided by the sum rules 

/ dx u v (x) — 2 and / dx d v (x) — 1 (4) 
Jo Jo 

to get the correct quantum numbers of the proton (and 
similarly for other hadrons). The gluon distribution is 
assumed to have the same basic Gaussian shape as the 
valence quarks, since they are all confined in the same 
region. The gluon normalisation is given by the momen- 
tum sum rule, i.e. 'YlnJndx xfi(x) = 1, where the sum 
also includes sea partons that are generated as follows. 

Sea partons arise from the non-perturbative dynamics 
of the bound state nucleon, for which it is appropriate to 
use a hadronic quantum mechanical basis. Therefore we 
consider hadronic fluctuations, e.g. for the proton 

\p) = a Q \p Q )+a pv \p7r )+a nlT \mr + ) + . . .+a AK \AK + )+. . . 

(5) 

Probing a parton i in a hadron H of such a fluctua- 
tion (Fig. \I]p) gives a sea parton with light-cone fraction 
x = xh Xi of the target proton, i.e. the sea distribu- 
tions are obtained from a convolution of the momentum 
K of the hadron and the momentum k of the parton in 
that hadron. The momentum K of the probed hadron 
is given by a similar Gaussian as in eq. J3J), with a sepa- 
rate width parameter ajy, and the momentum K' of the 
other hadron is then fixed by momentum conservation in 
the nucleon rest frame. Both hadrons are taken on-shell, 
which fixes their energies. This implies that energy is 
not conserved at this intermediate stage, but is of course 
restored for the observable final state. With the hadron 
four-vectors specified one obtains the light-cone fraction 
x H = K+/(K + K>) + . 

The above model for valence distributions is then ap- 
plied to give the parton momentum in H such that 
Xj = k + /K + is obtained. The flavour sum rules in 
eq. (0J must, of course, be modified to apply for H. 
The kinematical constraints to be applied in this case are 
m i < j 2 < xhW 2 and that the remnants (see Fig. [IJa) 
have invariant masses larger than their contained parton 
masses. A Monte Carlo method is used to simulate this 
two-step process by choosing K and k, impose the con- 
straints and obtain the momentum fraction x. By iterat- 
ing the procedure the sea quark and gluon distributions 
are generated. 
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The normalisation of the sea distributions is given by 
the amplitude coefficients a. These are partly given 
by Clebsch-Gordan coefficients, but depend primarily 
on non-perturbative dynamics that cannot be calculated 
from first principles in QCD and are, therefore, taken as 
free parameters. 

This model provides valence and sea parton x- 
distributions as shown in Fig. [21 These apply at a low 
scale Q 2 ,, an( i the distributions at higher Q 2 are obtained 
by applying the next-to-leading order QCD evolution 
equations as implemented in |l4| . The proton structure 
function F?(x, Q 2 ) can then be calculated and the model 
parameters be fitted to data from deep inelastic scatter- 
ing (DIS) resulting in the values: 

a u = 180 MeV, a d = 150 MeV, a g = 135 MeV 
<t h = 100 MeV, a 2 ea = 0.06, Ql = 0.6 GeV 2 (6) 

where a 2 ea is the fraction of the proton momentum car- 
ried by sea quarks at the scale Q\. This inclusive data 
can only be used to determine the overall normalisation 
a 2 Nv for the dominating light quark sea from fluctuations 
with pions in eq. (J5J . As shown in ^lj i the model repro- 
duces the Fi data well, which is remarkable in view of 
the model's simplicity with only six parameters. 

The model also reproduces the observed asymmetry 
between the u and d distributions as a result of the sup- 
pression of fluctuations with a ir~ relative to those with a 
7r + , since the former require a heavier baryon {e.g. A ++ ) 

For the NuTeV data the most interesting fluctuation 
is \AK), where an asymmetry arise since the s quark is 
in the heavier A baryon and the s is in the lighter K 
meson giving a harder momentum distribution for the s 
than the s. The same effect is present in fluctuations like 
|£AT) and \AK*), which are all implicitly included in our 
generic A A) fluctuation. Fluctuations where the ss pair 
is part of a meson wave function are here neglected since 
they are suppressed (<fi due to large mass and rj due to 
Clebsch-Gordan coefficients). 

The result of the model is shown in Fig. [5Jd, which 
clearly shows a difference between the xs{x) and xs(x) 
distributions as expected. Perturbative QCD evolution 
to larger Q 2 shifts these distributions to smaller x and 
gives rise to a symmetric ss sea arising from g — > ss, but 
the characteristic difference between xs(x) and xs(x) is 
still present as shown in Fig. [2t. 

The normalisation of the |AA) fluctuation cannot be 
safely calculated and is therefore taken as a free param- 
eter which was fitted to data from the CCFR collab- 
oration [l6| on the strange sea (s + s)/2, as shown in 
Fig- El The data can be reproduced with an ss contribu- 
tion from | A AT) fluctuations with normalisation such that 

Jq 1 dx(xs(x) + xs (x)) I Jq dx(xu(x) +xd(x)) « 0.5, i.e. the 
strange sea momentum fraction (at Qq) is approximately 
half of that of a light sea quark . in agreement with the 
parton density analyses in @, Il7| . This normalisation 
means that the coefficients a 2 in eq. JSJ scale essentially 
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Figure 2: Parton momentum distributions xfi(x) for different 
parton species i in the proton as obtained from the model with 
Gaussian momentum fluctuations in the proton giving valence 
distributions and in hadronic meson-baryon (e.g. AK + ) fluc- 
tuations of the proton giving sea partons. (a) All partons 
and (b) enlarged plot of the strange quark (s) and antiquark 
(s) sea at the scale Qo ~ 1 GeV 2 . (c) Strange quark and 
antiquark sea evolved with perturbative QCD to the scale 
Q 2 = 20 GeV 2 representative for the NuTeV data. Note the 
harder momentum distribution of the s quark (from A) rela- 
tive to the s (from K + ). 



with the fluctuation time At ~ 1/AE. Furthermore, 
this ss suppression is of similar magnitude as the ratio 
P(ss) / P(uu) k, 1/3 of probabilities for quark-antiquark 
production in phenomenological hadronization models, 
such as the Lund model [18|. Given the fact that both 
cases concern ss pair production in a non-perturbative 
process in a color field, this need not be surprising but 
indicate common features. 

With this strange sea we show in Fig. 0] the resulting 
asymmetry s~(x) — xs(x) — xs[x) and its combination 
with the folding function F(x) provided by NuTeV 
to account for their analysis and give the shift in the 
extracted value of sir? 0w- We obtain the integrated 
asymmetry S~ — f dxs~(x) = 0.00165, and the shift 
Asin 2 6W = J,, 1 dxs~{x)F{x) = -0.0017. Thus, the 
NuTeV value sin 2 9 W = 0.2277 ± 0.0016, which is about 
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Figure 3: CCFR deep inelastic scattering data |l6| on the 
strange sea distribution (xs(x) + xs(x))/2 in the nucleon at 
different Q 2 compared to our model based on \AK) fluctua- 
tions (with normalisation 's/u'= (s + s)/(u + d) as discussed 
in the text) at a low scale and evolved to larger Q 2 with per- 
turbative QCD that also adds a symmetric perturbative as 
component. 







— s(x) 
X£MF(x) 




L 





0.1 0.2 0.3 0.4 0.5 



Figure 4: The strange sea asymmetry a~ (x) = xs(x) — xs(x) 
(at Q 2 = 20 GeV 2 ) from the model and combined with 
the function F(x) accounting for NuTeV's analysis giving 
A sin 2 6 W = Jo dxs~{x)F{x) = -0.0017. 



3er above the Standard Model value 0.2227 ± 0.0004 Q, 
would be shifted to 0.2260 which is only 2. Oct above the 
Standard Model value. 



We note that the study in , based on exponentially 
suppressed AK fluctuations, gave an inconclusive result 
on the sign of S~ , and a too small strange sea to af- 
fect the NuTeV anomaly. In |2(| a result similar to ours 
was obtained. Unfortunately, these studies provide no 
comparisons with the measured strange sea, making the 
significance of their results difficult to assess. As shown 
recently |2lj higher order perturbative effects give a neg- 
ative contribution to S~ , although significantly smaller 
than our positive non-perturbative effect. 



The experimental situation is at present unclear. In 
0, a positive S~ was favoured based on several earlier 
experiments, whereas NuTeV 0, obtains the opposite 
sign based on their own data. However, the global anal- 
ysis in of the s — s asymmetry, including the NuTeV 
data as well as the CCFR data and using a very general 
functional form for s~(x), gives a best fit value for the 
asymmetry S~ of the same magnitude and sign as ours. 



We conclude that an asymmetry in the strange quark 
sea of the nucleon, which arises naturally in models where 
non-perturbative sea quark distributions originate from 
hadronic fluctuations of the nucleon, may reduce the 
NuTeV anomaly to a level which does not give a signifi- 
cant indication of physics beyond the Standard Model. 
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